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Introduction 


Sagan and Pollack (1967), Pollack and Sagan (1968), and 
Hansen and Cheyney (1968, 1969) have pointed out that 
some physical properties of clouds, such as mean particle 
size and optical depth, can be inferred from an analysis 
of the clouds' near infrared reflectivity. Belov/ we 
analyze observations of terrestrial water and ice clouds with 
this expectation in mind. It affords us an opportunity both 
to deduce the properties of these clouds and to check the 
consistency of our approach when applied to clouds we know 
something about . It is especially important to obtain 
tests of the plane-parallel approximation for real atmos- 
pheric cloud systems and the Mie scattering approximation 
for ice crystals. 

After outlining our computational scheme, we summarize 
the theoretical results for single scattering. These illus- 
trate the effects on the scattering of various input para- 
meters . Subsequently we discuss the infrared observations 
of Blau, et al. and attempt to remove the effects of gaseous 
absorption. In the next section we summarize the absorption 
corrections and relate them to the clouds ' temperature and 
to the mode of line formation. Finally, we estimate the 
phase , mean particle size , and optical depth of selected 
clouds from their spectral and angular scattering behavior. 



2 . Computational method 

To obtain the spectral reflectivity of clouds and the 
angular distribution of the scattered light we first com- 
puted the single scattering from a small volume containing a 
representative distribution of particle sizes and then the 
multiple scattering for the entire cloud. The single scattering 
computations were made using the Mie theory in essentially the 
way described by Deirmendjian and Clasen (1962) . However, the 
logarithmic derivatives of certain spherical Bessel functions 
occurring in the series were obtained by proceeding from the 
highest order terms? Kattawar and Plass (1967) have pointed out 
that an upward progression for the recursion relations is basically 
unstable and may cause significant errors for large size parameters 
The validity of our single scatter ing solutions was checked by 
making comparisons to published results of Deirmendjian (1964), 
unpublished company reports by the same author and unpublished 


work of H. Cheyney. 


We employed Deirmendjian 1 s (1964) "Cloud model" for the 
distribution of particle sizes, 

n(r) <x r exp (-6r/r m ) (1) 

where n(r) is the volume concentration of radius r and r 
v ' m 

is the radius at which the distribution has its maximum. 
Although sample size distributions from actual clouds often 
differ markedly from Deirmendjian’ s model, it serves the 
purpose of averaging out most of the large fluctuations which 
occur in the phase function (scattering diagram) for a single 
sphere . We have made preliminary computations to find the 
effect of changing the shape of the size distribution; the 
results for several different distributions indicate that the 
volume extinction . -{cr ^.) the single scattering albedo (o>) , 
the asymmetry factor (<' cos 9 )>) and the shape of the phase 
function (outside the region of the glory) depend mainly on 
the mean particle radius for extinction, 

? 

7 = 6 g ext < r ) n ( r ) rdr 

CD * ^ 

i °ext ( r ) n ( r ) dr 

ana not on the shape of the size distribution. This indi- 
cates that our computations in this paper with a single type 
of size distribution are meaningful, although in cases where 
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a more exact knowledge of the size distribution is available 

it should of course be used. For this paper our computations 

were made with r m - 2 , 4, 8, 16 and 32 ^ with the integrations 

over particle size extended to the radius r max = 25p for the 

first four distributions and to r = 50u for the case r ■ = 

max K m 

32|j . The integration increments were chosen small enough to 

make the phase function smooth (see, e.g., Dave, 1969a, 1969b). 

The values used for r were somewhat arbitrary but that is 

max r. 

not essential since the purpose of the integration over particle 

sizes was to smooth out single particle effects; we usually 

found r’ — r . 

m 

Except where otherwise indicated we employed the optical 
constants for water and ice tabulated by Irvine and Pollack 
(1968) ; the computations were made at each wavelength (~ 50) at 
which the authors tabulate the optical constants for the region 
of interest (1.2 > \ > 3.6[j.) . 

The intensity of light multiply scattered by the clouds 
was obtained by using the "double only" computing method, 
described by Hansen (1969a), which is a variation of van de 
Hulst's (1963) doubling method. If errors ^ 1% are tolerable, 
as is certainly the case for comparison to most observations, 
then the computing time may be greatly speeded up. Some of 
the more useful shortcuts which we have tested are described 
in the appendix. ’ -/ ; ~v 
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3 . Results for single scattering 

The single scattering behavior of an ensemble of aerosols 
is a function of the complex index of refraction, n c “ n r “ *- n i 
and the distribution of particle sizes, expressed in units of 
the ratio of the Circumference of the particle to the wave- 
length x = 2nr/\ . Assuming that the distribution of size 
is given by (1) , we find that the parameter x m , which equals 
2 TT r m /^ > serves to define the dependence of the scattering 
behavior upon particle size. The results described below are 
of both general interest and of help in understanding the 
infrared properties of terrestrial clouds analyzed in later 
sections. In the calculations for Figs. 1-3 the integrations 
over the size parameter x were performed to an upper limit of 

x = 2x , except when x = 128 for which x was set equal 
max m m max 

to 200 . 

a. General single scattering results 

Fig. 1 shows the dependence of the normalized phase func- 
tion p(0) (Deirmendjian, 1964, 1969) upon the particle size 

parameter x in the case of no absorption (n. =0) with the 
m r 

real refractive index equal to 1.33, the value for liquid water 


in the visible. The variable 0 is the angle of scatter. 

The individual curves have been vertically displaced 
from one another. The short horizontal line intersecting 
a given phase curve denotes the position at which the phase 
function has a value of unity. For the largest values of 
x m , the very precipitous decline in the value of p near 0° 
corresponds to the diffraction peak (Deirmendjian, 1964, 

Dave, 1969c) ; the strong maximum at a scattering angle of v 
142° is the main rainbow (caused by rays undergoing a single 
internal reflection) with its first supernumery bow (van de 
Hulst, 1957, p . 241) located at 147°; the second rainbow 
(2 internal reflections) is located at 123° while its first 
supernumery bow lies at 114°; and finally the glory corresponds 
to the overall increase and oscillatory behavior near 180°. 

As x m decreases all these features become less pronounced and 
broader; in addition the rainbow shifts in location towards 
larger scattering angles and the? slope of the phase function 
decreases until near x m = 1/2 it is very similar to the Rayleigh 
phase function. 


The effect of introducing absorption within the particles 
is investigated in Pig. 2. For x m = 32, the glory and rain- 
bow have been effectively suppressed when n^ is comparable 
to or larger than 0 . 03 . This effect can be understood using 
concepts from geometrical optics. A ray traversing a path 
equal to the particle radius will be diminished in strength 
by exp(-ka) = exp (-2 x n^) where k is the linear absorption 
coefficient . Setting n^ equal to 0 .03 and x to 32, we see 
that the ray is diminished by almost a factor of 10 in in- 
tensity. Similarly the diffraction component which is not 
affected by absorption, extends to larger scattering angles 
as nj_ increases. 

For x m again of 32 , the value of the phase function in 
the backward hemisphere (9 > 9 °°) begins to increase as n. 

JU 

becomes larger than 0 .1 . This may be attributed to an en- 
hanced value of externally reflected light, as seen from the 
Fresnel equation. For the smaller particles (x m =2) the 
results are basically similar, although less pronounced . 

Finally we study the influence of the real part of the 
index of refraction in Fig. 3 . For x m = 32, increasing n r 
causes the rainbow to shift to larger scattering angles, 
until it merges with the glory, resulting in a large increase 
in p (l80°) . Aside from the diffraction peak, small angle 
scattering becomes more dominant as the refractive index 
decreases, a result which shows up in the asymmetry factor 
< cos 0 > described below. We will find this effect to 
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be of significance in understanding certain spectral charac- 
teristics of water clouds. Similarly at the lower values of 
n r refraction tends to dominate over diffraction at smaller 
angles of incidence and hence the break in the diffraction 
peak occurs at smaller values of 0 . The right hand half of the 
figure shows that similar effects occur for particles with 
absorption (n^ = 0.01) and small particles (x m - 2) . 
b. Single scattering for water and ice 

One important integral scattering parameter is the single 
scattering albedo, uJ, the ratio of the amount of light scattered 
to that which is scattered and absorbed . In Figs . 4 and 5 the 
spectral variation of 1 - uj is exhibited for liquid water and ice 

particles, respectively, for five values of the size’ parameter r m . 

1 

' ' ' ' ....... . ^ 

As noted by Irvine and Pollack (1968) , the maxima and minima 
of ice and water are displaced by about 0 .lp, in wavelength , a 
feature useful in distinguishing the phase of water clouds. Also 
the single scattering albedo generally declines monotonically with 
increasing particle size. It should be noted that in all cases 
x m > 3 . For sizes much smaller than this the particles will become 
completely absorbing and the above generalization is no longer true 
(van de Hulst, 1957) . An interesting mild deviation from the 
general variation of uu with particle size occurs near a wavelength 
of 3ji and may be due to the perturbing effect of the smallest 
particles having values of x < 1. Near 3p the value of n^ is so 
large that the particles become completely opaque so that allowing 


for the diffraction peak & is about 1/2. 

The asymmetry factor < cos 9 > is defined an the solid 
angle average of cos 9 weighted by the phase function. It 
describes the degree of forward scattering of the phase 
function. For isotropic scattering it has a value of sere 
while it approaches unity as small angle scattering tends to 
dominate. Figs . 6 and 7 show the spectral behavior of 
< cos 9 > for liquid water and ice particles, respectively. 
For the largest value of x m there is a very pronounced 
increase in <( cos 9 )> slightly shortward of 7a and a less 
obvious minimum somewhat longward of . This behavior 
reflects the anamolous de spersion change s in n near the 
very strong absorption feature centered near pp, with changes 
^ also influencing < cos e > near its miniwa. As shown 
in Fig. 3 and discussed earlier the phase function becomes 
more forward scattering as n r approaches unity. 

From multiple scattering computations reported below 
and elsewhere (Hansen, 1969a; Hansen and Cheyney, 1968) 
we find that for wavelengths less than 2.gp and between 3.-4 
and 3.6a the reflectivity of thick clouds depends primarily 
on the single scattering albedo 1 qualitatively the curves 
of log (1 - 35) are very similar to the spectral variation 
of the cloud reflectivities . As the single scattering albedo 
decreases, the .less probable it is for a photon to survive 
a number of scattering events. Since the single scattering 


10 


albedo in these wavelength regions varies systematically 
with the characteristic particle size r m (cf. Figs. 4 and 5) 
we can obtain particle size information from the near infra- 
red spectral behavior of clouds, as has been pointed out by 
Sagan and Pollack (1967) , Hansen and Cheyney (1968) , and 
Irvite and Pollack (1968) . 

On the other hand in the spectral region between 2.3 
and 3.4a the single scattering albedo has a constant value 
of about 1/2 and the spectral behavior of the cloud re- 
flectivity depends primarily on <( cos 9 )>. There is an 
inverse qualitative similarity between the wavelength varia- 
tion of <( cos © ]> and computed cloud reflectivities. As 
<( cos 9 y increases , less light incident on a cloud layer is 
reflected back out in the first few scattering events . If 
the single scattering albedo is sufficiently Je ss than one 
so that these scattering events are the major producers of 
reflected photons , the cloud reflectivity will decrease as 
<( cos 9 )> increases. For example the maximum in < cos 0 > 

slightly shortward of 3q results in a minimum value for the 
cloud reflectivity at the same wavelength, as shown below. 
The spectral variation of < cos 0 > is a reflection of the 
spectral behavior of the indices of refraction. 


4 . Reflectivity measurements of terrestrial clouds 

i 

Blau, Bspinola, and Reifenstein (1966; see Blau and 
Espinola, 1965 » for a more detailed report) have obtained 
airborne infrared measurements of the reflectance properties 
of terrestrial clouds, These consist of both spectral and 
angular observations. Observations were performed from an 
aircraft above the cloud of interest and on a given day- 
measurements were made either in the 1 .2|i to 2 , 5 \x region or 
the 2.4jJt to ?. 6 \i wavelength domain. In addition on some 
occasions angular scattering information was obtained by 
flying along a hexagonal path and performing spectral measure- 
ments of the cloud area situated at the center of the hexa- 
gon. In such measurements the angles of reflection and 
incidence remain constant, while the azimuth and angle of 
scatter vary. Below we describe these observations in greater 
detail and outline the transformation we applied to them so 
as to be able to compare them, with our multiple scatlerlhg^^^^^^^^^^^^^^ 
computations . 

The spectral radiance (specific intensity) values 
reported by Blau et al. refer to averages of a number of 
spectra obtained close together in time. In addition they 
also give values for the standard deviation of each average 
value . The standard deviation is not necessarily a reflec- 
tion of the error of measurement, because it also includes 
short term variations in the properties of the observed cloud . 
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¥e subjected these measurements to two types of transforma- 
tions . First we divided the observed specific intensity by 
the solar flux outside the atmosphere of the earth, F^: 

R * = Vtt “ 1C FT (5) 

where S, defined by Eq. (3) , is the usual scattering function, 0 O 

cos 1 |j 0 the angle of incidence, and 0 = cos” 1 ^ the angle of 

reflection. If the clouds were a lambert surface, i.e., 

if they scattered Isotropically, the normalized reflectivity 

would equal the spherical albedo of the clouds and would 

exhibit no angular variation. 

All the spectra show absorption features due to water 

vapor. In addition some absorption due to carbon dioxide 

at 2p is expected and there are strong COg absorption features 

near 2 »7p. Because gaseous absorption takes place not only 

above the clouds but through multiple reflection within the 

clouds, we did not exactly correct for this effect. One 

cannot practically correct for the multiple scattering within 

the clouds by assigning an effective single scattering albedo 

to the gaseous absorption component, because the spectral 

resolution was much larger than the spectral domain over 

which the gaseous absorption is constant. We performed our 

water vapor absorption corrections by comparing the value of 

R at two spectral positions, one position expected to have 
A 
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very little gaseous absorption and the other a large amount. 
Furthermore the clouds were expected to have approximately 
the same intrinsic value for at the two wavelength posi- 
tions. In the 1.2 - 2.5p region we compared values at 
1.28 ^i with either 1 . 35 p or 1.58 p, while for the 2.4 - 3.6|x 
domain we considered values either at 2.50 or 2.55 with ones 
at either 2.60 or 2 . 6 lp. The ratio of the values of R^ were 
equated to ratios of gaseous transmissivities calculated by 
Wyatt, Stull and Plass (1962) . The comparison was made at 
pressures and temperatures closest to the cloud top condi- 
tions, as inferred from the cloud top altitude, and in all 
cases at an effective re solution of 100 cm” ■*". While this 
resolution is somewhat poorer than the resolution of the 
spectrometer it was found to yield the most consistent 
results . In part the need to employ lOOcrn”^ resolution 
is a reflection of the breakdown of the theoretical statistical 
model at finer resolutions, as indicated by the appearance 
of high frequency features whose amplitude is too large . 

From the comparison with Wyatt et al. ' s tables we derived 
an effective water vapor abundance , W, and used this amount 
to estimate the effective gaseous transmission at other 
wavelengths . At wavelengths where there is more absorption 
by the cloud aerosols , there is less multiple scattering and 
in this sense the transmission correction is an overestimate . 
This circumstance pertains at most of the other wavelengths 



and the true cloud reflectivity in general should lie between 
the uncorrected and "corrected" values. 

Absorption effects by carbon dioxide were much more 
localized in the spectrum. For the strong absorption band 
near 2 . 7 M-* we derived an effective amount of CO^ and an 
average pressure by allowing for the path down to and up from 
the cloud top and estimating the effective path in the clouds 
from the water vapor absorption amount. The latter can be 
related to the atmospheric temperature where the absorption 
takes place, as discussed below, and to the pathlength with 
the help of the U.S. Standard Atmosphere Supplements (I966) . 
For the much smaller corrections within the 2 \x band, we used 
the water vapor absorption amounts as a guide in a less rigor- 
ous fashion. Transmission corrections were then obtained 


from the tables of Stull, Wyatt and Plass (1965) . 

Fortunately the angular measurements refer to the path- 
lengths of nearly constant angle of incidence and reflection 
and so to first order they require no absorption correction. 
However because some absorption takes place through multiple 
scattering the effective atmosphere transmission may have an 
azimuthal dependence. This effect is very difficult to correct 
for and no attempt was made to do so. 


***** 


5. Discussion of gaseous absorption corrections 


In Table 1 we summarize the water vapor amounts, W eff , 
deduced from Blau et al. f s spectra. The first column gives 
the figure number of the spectra in their final report. 

P is the pressure assumed in deriving W Also given are 

the cloud type and the cloud top altitude. Using the U.S. 
Standard Atmosphere Supplements ( 1966 ) , we have estimated 
the temperature at the cloudtop, T c( p from its altitude. 

To assess the contribution of multiple scattering within 
the clouds to the observed water vapor absorption features, 
we have computed the equivalent amount of water vapor W 
which the sunlight passes through above the clouds, on its 
path down to the cloudtop and up to the airplane. In per- 
forming the calculation we have assumed the atmosphere to be 
saturated . In addition we have corrected for the difference 
in the value between the pressure used in obtaining 
and the actual cloud top pressure by assuming that pressure 
and gas amount are equally effective in causing absorption. 
For the first three spectra of the 1.2 - 2.5u region W ^ 
is significantly smaller than W . Hence most of the absorp- 
tion takes place in the clouds. A confirmation of this de- 
duction is obtained by comparing the values of W eff for the 
two sets of spectra from Fig. 55. The value of obtained 

for the larger angle of incidence is smaller than that for the 


TABLE 








smaller angle of incidence; this finding is opposite to the 
expectation for absorption taking place above the cloud, but 
in accord with predictions for multiple scattering within a 
cloud layer (Chamberlain, 1965) , Such a phenomenon is also 
present for at least some of the gaseous absorption features 


of Venus (Chamberlain and Kuiper, 1956). 

On the other hand, W ref is larger or comparable to 
for the two spectra pertaining to the 2,4 - 3.6p region. 


This results in part from the eloudtops being located high 


within the stratosphere where the water vapor abundance is 
only a few percent of the saturation abundance. For these 


spectra we can conclude that the fraction of the cloud signi- 
ficantly contributing to the scattering lies within the 
stratosphere and therefore we are surely viewing ice 
particles . The relatively small depth of penetration for 
these spectra is a result of the highly absorbing nature 
of ice aerosols at 2 and 2 . 6q , the wavelengths at which 


W eff was obtained. 

To derive an estimate of the depth of penetration for 
the other clouds, we have computed T ^ , the base tempera- 
ture required within a saturated atmosphere so that light 
traveling on a straight line down to and up from this level 
at a 60° angle would experience the observed amount of absorp 
tion. Since the actual pathlength within the clouds is more 


tortuous 5 T , is probably a slight overestimate of the 
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level of line formation. We see that for the first three 
spectra T g ^ is substantially larger than T d , implying 
substantial penetration within the clouds. At the wave- 
lengths used to derive W eff> for these spectra the cloud 
aerosols are essentially transparent. 


6. Analysis of the cloud spectra 

Our theoretical spectra are functions of three parameters ; 
the characteristic particle size r. the optical thickness t 
of the clouds, and the phase of the cloud, i .e . , whether the 
aerosols are liquid water or ice* Below we attempt to 
estimate each of these parameters by comparing the theoreti- 
cal and observed spectra . 

In Figs . 8 and 9 Blau et al's observations of a cirrus 
cloud in the 1.2 to 2.5^ wavelength region are compared with 
theoretical spectra for an ice cloud. The observations 
correspond to Blau et al . ’ s Fig . 55 for an angle of incidence 
of 71°. Circles and solid bars represent the average values 
and standard deviations of the reflectivity R (eq. 5) after 
correction for gas absorption, while the uncorrected observa- 
tions are indicated by triangles and dotted bars in the cases 
where they differ significantly from the corrected values . 

As explained above, the absorption corrections may be over- 
estimates in regions where the cloud aerosols strongly absorb. 
In this case the true reflectivity will be somewhere between 


the corrected and uncorrected values. The theoretical curves 

in each portion of the diagrams correspond to various choices 

for the optical depth; the particle size parameter is vailed 

by a factor of 2 between successive sections . Comparing the 

four sections of Figs. 8 and 9, we see that the observations 

permit the determination of a well defined value for r m , l6ji. 

In all of the figures the optical depth refers to A = 1.2|i 

however, the wavelength dependence of t is small since we 

always consider a distribution of particle sizes with x m >1 . 

In determining t, values of R between 1.2 and 1 . V are of 

particular use : the theoretical curves differ at these 

wavelengths by a maximum amount . For- the cloud observations 

illustrated in Fig. 9 a value r ~ 10 appears to give the 

best fit . However , the optical depth is more difficult to 

estimate than the particle size, and this derived optical 

depth should probably only be regarded as a lower limit. 

Danielson, Moore and van de Hulst ( 19 69 ) hare argued , from 

observations and computations, that condensation nuclei 

limit the cloud reflectivity at wavelengths where ice and 

-■3 

water do not effectively absorb [1 - 0 <,10 .] . 

As mentioned above the absorption features of water and 
ice are displaced somewhat from one another. For example 
peak absorptions and hence minimum values of R occur at 1.45 
and 1.95p for water and at 1.52 and 2 .OOp, for ice . As a 
result the theoretical water spectra do not fit the observed 


spectra in these regions. On the other hand, there is some 
indication that a local maximum in R occurs at 2.2(j, a posi- 
tion expected for water clouds, rather than 2 .3|j, the place 
for ice clouds. A similar situation was found in some of the 
other spectra. This could be understood in terms of a mixed 
phase model with ice predominating near the top of the cloud, 
but the quality of the present data did not warrant calculations 
for such a model . 

In a similar fashion we analyzed several other observed 
spectra. The results are summarized in Table 2. It is en- 
couraging to see that the same cloud properties were found from 
observations of the same cloud viewed at two different solar 
illumination angles (Fig. 55) . An interesting feature of 
Table 2 is the similar value for r m found for the various ice 
cloud spectra. The deduced average particle radius of 16|i, or 
diameter of 32^, is consistent with values typically obtained 
from direct sampling measurements of ice clouds. 


TABLE 2 

DEDUCED PROPERTIES OP SOME OBSERVED TERRESTRIAL CLOUDS 


Bl&u * s 



Deduced 

Pror 

)erties 

Pig. 

9 0 

Cloud 

Phase 


1 T 



X * 1.2 to 2.5u 




55 

71 

cirrus 

ice> 

16 

10 

55 

64 

cirrus 

ice 

16 

10 

55 

.44 

cumulus 

water 

8 

4 

51 

52 

thick cirrus 
over cumulus 

ice 

16 

10 



A = 2.4 to 3* 6u 




18 

48 

cumulus 

ice 

16 

>5 

20 

52 

cirroform 

ice 

16 

>5 


In Pig. 10 we consider the degree to which multiple 
scattering occurs in the 2.5 to 3.5li spectral region, 
where the cloud aerosols are highly absorbing. The curves 
illustrated are for an optical depth of unity and a large 
enough optical depth (t = 128 ) to be considered equivalent to 
an infinite value; we obtained almost identical results for 
all optical depths in excess of 10 . The dashed curve repre- 
sents photons scattered only once , the dotted curve photons 
scattered n times (n > 1) with the first n-1 scattering events 
occuring within the diffraction peak (Hansen, 1969b), while 
the solid curve represents all the photons reflected from 
the cloud. We see that even if the dotted curve is considered 


as representing single scattering, single scattering computations are 
inadequate to describe the reflectivity in this highly opaque part of 
the spectrum. 

In Fig. 11 we compare computed cloud spectra for t s 10 with data 
given in Blau et al.'s Figs. 18 and 20. The theoretical curves corre- 
spond to various choices of r m . The data points between 2.5 and 2 . 6(j 
indicate a particle size of about 16 jj in both cases. Near 3jj the ob- 
servations and calculations are not in good agreement; this disagreement 
is discussed in detail, in $ 8. 

Finally we note that sample calculations of blackbody thermal emis- 
sion from the clouds showed this contribution to be very small compared 
to reflected sunlight, even at the longest wavelengths of observation. 

7. Analysis of the clouds' angular scattering 

In this section we attempt to assess the information content of 
the clouds ' angular scattering properties. As mentioned earlier, ob- 
servations were made at nearly constant angles of incidence and re- 
flection but with varying azimuth and hence scattering angles . In 
Figs. 12 and 13 we contrast the theoretical single scattering behavior 
of water clouds with the complete multiple scattering behavior. The 
computations were made for angles of incidence and reflection of 60° 
and 80°, respectively. In Fig. 12 the reflectivity has been calculated 
for a distribution of water particles having a mean size r m of 2^, while 
Fig. 13 pertains to a value of 16(j for r^. For the larger size particles 
we defined single scattering to include photons scattered (n-1) times 


in the diffraction peak (n > 1) before suffering a final 
scattering back out of the cloud . Since the 2\i size 
particles do not have a striking diffraction peak, we 
defined single scattering for them in the conventional 
manner. 

We see that the effect of multiple scattering is to wash 
out features such as the rainbow and to greatly diminish the 
angular variation of R. Even for an optical depth of unity 
the actual scattering behavior is markedly different from 
that of single scattering. At an optical depth of 16 the 
scattering behavior is quite close to that for an infinite 
optical depth. On the basis of their high albedo in the 
visible we know most clouds have an optical depth of 16 or 
more . In this event a knowledge of the exact value of the 
optical depth is not too important for being able to predict 
the reflectivity properties. 

In Fig. 14 we compare theoretical calculations with the 
angular dependence exhibited by a cumulus cloud whose top 
was at 4 kft . This doud was part of Hurricane Gladys. For 
each data point the angle of reflection was 80° but the 
angle of incidence varied between 56° and 64°. The 
theoretical computations were performed for spherical water 
particles at the angles appropriate for each observation, 
and the theoretical points were joined together by straight 
lines. As mentioned above, no correction for gas absorption 


was made . 
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All the theoretical .spectra are normalized to fit the 
data point at the lowest value of (0 - CL) . A particle size 
of 8{jl appears to yield the best fit to the measurements, a 
result compatible with typical direct sampling measurements 
of cumulus clouds . Unfortunately, the absorption corrections 
were so large we were unable to meaningfully analyze the 
spectra to confirm this deduction. 

Angular computations in the 2 .4 to 3.6p region are shown 
in Figs. 15 and 16 , The calculations were made for 
spherical ice particles with t he realization that this may 
lead to a very poor approximation to the true angular scatter 
ing behavior of ice clouds . In Fig. 15 a mean particle 
radius r ■ ' of l6|_t wasused and the optical depth varied. 

In Fig. 16 the particle size is varied, while the optical 
depth was set equal to 32. In all cases the angles of 
incidence and reflection are 50° and 80° respectively. 

We see that for cases of intermediate aerosol absorption, 
e .g. , r m - l6|i and % = 2.72p, the exact solutions most 
clearly preserve such features as the rainbow peak and 
exhibit the strongest dependence upon particle size. When 
the aerosol absorption is relatively small, many scattering 
events occur leading to a smooth, featureless angular 
behavior. When the aerosol absorption is very large the 
rainbow is not present even for single scattering. 

In Fig. 17 we compare theoretical and observed cloud 
reflectivities for an angle of reflection of 80° and an 
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angle of incidence between 44° and 56°. The observations pertain to cir- 
roform clouds whose tops were at 43 kft. They were part of Hurricane 

I . ■ . . ' • . • • • •• • ’ " • • • ■ ■ 

t . • . • •' • - .••••••••' ' . ' •• • . ' 

Gladys . The theoretical curves have the same meaning as those in Fig. 14, 
except that Fig. 17 is for spherical ice particles. A particle radius of 

i 

16|i, a value found typical for cirrus clouds from the spectral analysis, / 

gives no worse a fit to the data than other sizes. The calculations at 

all 3 wavelengths were made with the optical constants given by Irvine 

and Pollack? however, if we accept the modification in the optical constants 

which is indicated by the results in § 8, then the theoretical calculations 

at \ = 3 .10|i would be changed and brought into better agreement with the 

observations. At 3.10(j the primary effect of the new optical constants 

is to increase the Fresnel reflection in the backward direction? this 

would increase the theoretical reflectivity at the two right-most data 

points in Fig. 17. 

8. Optical constants of ice near \ = 3|j. 

In this section we reconsider the ice cloud spectra near X — 3 p where 
the calculations and observations are not in good agreement, and we sug- 
gest one possible explanation for the discrepancy. i 

A minimum in the observed spectra (Fig. 11) occurs at about 2.9p in 

wavelength and a maximum around 3.1p. These features are particularly 

• ■ ' ■ - . - | 

prominent in the right hand graph. This may at first appear surprising 

. . . \ ■. y.y y 7 . / : -v; yy.v • • y- • . i ; ' . . 

since ice has its maximum absorption at 3. Ip (Irvine and Pollack, 1968) * 

Blau and Espino la (1969) first pointed out the minima in the spectra and 
correctly attributed it to the anomalous dispersion of n r , which is dis- 

v 7".7 ■ ■■ ■ -.. ''7 k 7 ^ 7777 ^ ... ;/ yk7y;kk7 7.7 7^^^ 

cussed above. A local maxima in the reflectivity is expected near 3.1p 

for similar reasons. Our computed spectra qualitatively show these effects, j 

but particularly for the right hand graph (Fig. 11) they fail to quant it a- 


! - / . • •; .y: ••••. . y-..7-'- "v . * / . 'y y •'/ ' . ^ y ; ...y. ’ ' ' ' • . \ 

tively match the observations. Blau and Bspinola encountered a similar 
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difficulty in explainina the minimum and suggested that Irvine and 


Pollack's value of n r be revised near this wavelength position. 

We have also ccnsidered the possibility of revising 
Irvine and Pollack's optical constaits for ice in the vicinity 
of the strong 3. Ip. feature. We have attempted to do this in 
a consistent fashion by 'relating changes in the real part of 
the index of refraction, n r , to changes in the imaginary 
part n.j_» Spitzer and Kleinman (1961) have given relation- 
ships between n r and n^ under the assumption that individual 
absorption bands may be considered as classical oscillators 
and have obtained a very good fit to reflection measurements 
of strong infrared bands of quartz. Assuming that the 3 . Ip 
band is a dominant feature at nearby wavelengths and neglect- 
ing small differences between the value of the frequency at 
the desired positions and the central frequency of the band, 
we have simplified Spitzer and Kleinman' s formulas to the 
following general form: 

pF^ 

n r n l = (AA) 2 +r a W 



a + 

(AA) 2 +r 2 



where p Is the maximum value of the product of n n^, r one 
half the value of the width between half maximum points of 


n r n i' ^ the difference in wavelength between the value of 
interest and the position of the maximum, and A is a constant. 
These formulae give a good fit to the water and ice data 
of Irvine and Pollack near the 3u absorption feature. 


For the ice constants given by these authors p is about 0.4-2, 
and F about 0 . 15p . The constant A is readily found by apply- 
ing the second formula to wavelengths significantly short- 
ward of 3p; we find A is 1.66. 

In modifying the indices of refraction given by Irvine 
and Pollack we assumed that r was unchanged and varied the 
only remaining free parameter p by various scale factors c . 

We have modified the published data between 2. 8p and 5.6(1 
because this is the region dominated by the 3. Ip band . In 
addition the published indices shortward of 2 ,8p were based 
on a different set of measurements than those used at 2.8p 
and longward. These first set were checked against other 
measurements and found to be in good agreement . The results 
of two trial modifications are given in Fig. 18 along with 
the original values . In addition the effects of the modified 
indices on the single scattering albedo cs and the angular 
asymmetry parameter <( cos 9 y are shown. At short wavelengths 
the value c = 2 led to values of n r less than 1 and so this 
portion of the curve has not been drawn. 

Since the optical parameters have not been changed 
shortward of 2 .8p, the deduced characteristic particle sizes 
and the lower limit on the optical depth will still hold . 

We were able to deduce these from the data points between 
2.47p and 2.65|i. Fig. 19 shows theoretical spectra for Blau's 
Figs. 18 and 20 based on values of c = 1 (the old optical 
constants ) , 1.6, and 2. We see that a value of c of about 


1.6 leads to a good fit of the spectra for both Figs. 18 and 20. 

Some justification for our method of varying the ice constants 
is provided by the followings For wavelengths < 2.8ji the optical 
constants of Irvine and Pollack were based on several sources of 
data which agreed with one another. However, for \ ^ 2.8ja Irvine 
and Pollack had only one source for the ice absorption coefficients 
and this source gave water absorption coefficients for \ ~ 3|j which 
were ~ 60% less than the values obtained by several other experi- 
menters. Thus the proposed revision of the indices of refraction 
seems quite reasonable. 

9 . Summary 

This study indicates that characteristic broad band absorption 
features in light scattered by clouds can be used to help identify 
the scattering material and to determine the particle size and cloud 
optical depth. The results are consistent with the assumption that 
it is adequate to employ the spherical particle approximation in 
calculating the spectral reflectivity. As far as the angular be- 
havior is concerned, not enough accurate observations are available 
for an adequate analysis of the theory. More extensive measurements 
on atmospheric clouds are desirable, as well as carefully controlled 
laboratory observations along the lines of the recent measurements 
by Zander (1968) and Plummer (1969) . 
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APPENDIX 

The following shortcuts in the numerical work have been 
tested and found useful for the multiple scattering computa- 
tions : 

(a) For large particles (x = 2tt a/7v > 25) there is a 
sharp diffraction peak in the phase function which nec- 
essitates a large number of terms in the cosine expansion 
of the scattering function and a large number of points in 
the integrations over \x = cos © . However, in the case of 
conservative scattering we have shown elsewhere (Hansen, 1969b j 
see also Potter, 1969) that the photons scattered into the 
forward spike may be approximated as being unscattered by 
truncating the forward peak from the phase function and 
reducing the interaction optical thickness t such that 

; ; : 

where t’ is the optical thickness to be used with the trun- 
cated phase function, F is the fraction of photons scattered 
into the forward peak, 

s' = 5 (p - p’ ) (A2) 

4 7T 

and p and p’ are the untruncated and truncated phase func- 
tions, respectively. For conservative scattering this apprcximat 


Introduces large errors in the reflected intensities only 
for small total scattering angles (0 ~ 0°); it introduces 
errors of a few percent if the total scattering angle corre- 
sponds to a sharp feature in the phase function (such as the 
glory) and it also introduces errors of a few per cent if the 
incident or emergent angle is near grazing (9, 9 q , ~ 90°) 
elsewhere the error is 1$. 

In the case of nonconservative scattering the single 
scattering albedo* co, must be scaled such that 


1 



because the assumption that photons in the forward peak 
should be treat ed as unscattered implies tuat the absorption 
cross section is unchanged and the scattering cross section 
reduced by the factor 1 - F. Several tests of this approxima 
tion were made and it was found to become increasingly 
accurate as cd decreased. Even in the wavelength region ~ 3u 
where often ~ 80$ of the single scattered photons were cut 
off with the truncating of the diffraction peak, the error 
introduced was Vfo except for scattering angles ~ 0°. 

(b) In the expansion of the scattering function in 
cosines , we write: 

od 

S(t; p., p; u. , p ) = E S ra (T;.|i, M- 0 ) cos m (P - P a ) (A4) 

m=0 




52 

where U Q and U denote the cosine of the angles of incidence 
and reflection, respectively and 0 and 0 are the correspond- 
ing azimuth angles. The number of terms, M, needed to obtain 
an accuracy within 1 % for all \i and \x Q with a strongly 
anisotropic phase function is typically about one hundred. 
However, the range of u and u 0 for which the numerical value 
of S m (u, u 0 ) is not negligible decreases as m increases 
until for m ~ M only S(tj u ~ 0 \i Q ~ 0) is significant. 

If S (t; u, 0, u 0 > 0 O ) is to be calculated at N values of 

U on the interval (0,1) and N values of u Q then S 0 (tj u, U G ) 

2 

must be computed at N points and in the integrations over \x 
which include S°(t; \i, u 0 ) as a factor, N points are employed 
however the number of points required decreases steadily as 
m increases until for S M (t$ [l, u 0 ) calculations are only 
needed for one point and only one point is needed in the 
integrations . For a given accuracy specification it is 
easy by numerical testing to find the number of points at 
which S m (Tj \x, u c ) Must be calculated and included in the 
integrations , A factor of 2 - 3 in computer time may be 
saved with an introduction of errors ~ 1 - 2 %. 

(c) The strongest azimuthal dependence arises from 
single scattered photons but an analytic expression exists 
for the intensity due to these photons (Hansen, 1969a) and 
hence considerable computing time may be saved by writing 
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S(TJH,0Jii o ,0 o ) = S ss (TJ|i,0; 1 i oJ 0 o ) 

CD 

+ L CS m (Ti|x,jx ) - Sg (T;p.,^ ) ] cos m(0 - 0 ) 
m=0 



where the subscript ss labels the contributions of single 
scattering to the S function. Typically the number of terms 
needed in (A5) is ~ 500 of the number required with (A4>) 
to achieve the same accuracy. 

(d) Several additional ways to save computer time, 
which we found by numerical experimentation, can be shown 
to have a firm theoretical basis from work of van de Hulst 

( forthcoming book) . van de Hulst shows that each term 

in the cos m (0 - 0 Q ) expansion may be thought of as having 
an effective albedo for single scattering and this albedo 
decreases steadily as m increases . Some consequences are : 
the doubling process for terms with m > 1 may be initiated at 
an optical thickness t *\> 2’^ rather than 2"^- ? j for m > 1 

the asymptotic value of the scattering function is obtained 
already at t ~ 8; for m >M/4 the sum of the infinite series 
occurring in the doubling equations may be replaced by the 
value of the first term. These simplifications may easily 
reduce the computing time by a factor ~ 5. 

If all four of the above methods for reducing the com- 
puter time are employed the total time Caving is not the 
product of the factors which each gives alone because there 
is considerable overlap. In the computations for this paper 
(performed on an IBM 360/95 ) we always employed (d) and for 
the particle size distributions with r = l6p and J>2\i we 




i 3* 

i : 

I 

J. 

i used (a) . All of the above methods may be worthwhile for 

| \ 

! slower computers and especially for problems such as line 

formation and the multiple scattering of polarized light. 

Note added in revision: Dave (1970) has independently shown.. 

/ 

that the number of terms required in the Fourier expansion 

•i. ■ , .. . 

depends strongly on ^ and [shortcut (b) above] and he has 

presented graphical illustrations of this. Dave also makes use 

of the fact that the effective albedo decreases toward higher 

.... 

terms in the Fourier expansion [shortcut (d) ] 
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FIGURES 


Single scattering phase functions for a size dis- 
tribution of transparent spherical particles with 
a real refractive index typical of water and ice 
in the near infrared $ the curves show the effect 
of changing the characteristic particle size. 

In Figs . 1-5 the vertical scales apply to the 
uppermost curve on the left side and the scales 
for the other curves may be obtained by multipli- 
cation by a power of 10 such that the horizontal 
bar on each curve occurs at p(9) = 1. 

Single scattering phase functions for a size 
distribution of spherical particles showing the 
effect of absorption within the particles for 
large particles (x m = 32) and particles of 
moderate size (x m = 2) 

Single scattering phase functions for a size dis- 
tribution of spherical particles showing the effect 
of changing the real part of the refractive index 
for large particles with no absorption (left), 
for particles of moderate size with no absorption 
(upper right), and for large particles with moderate 
absorption (lower right) . 


Fig. 4 


Single scattering albedo for a cloud of spherical 
water particles for five different particle size 
distributions . 

Fig. 5 Same as Fig. 4 for spherical ice particles. 

Fig. 6 Asymmetry factor of the phase function for a cloud 
of spherical water particles for five different 
particle size distributions . 

Fig. 7 Same as Fig. 6 for spherical ice particles. 

Fig. 8 Theoretical cloud reflectivities for 9=0, 

0 Q = 71° and ft - ft 0 - 180° for five cloud optical 
thi ckne sses and two s iz e dis tribution s of 
spherical ice particles . The circles and solid 
bars represent observations by Blau, et al . of 
cirrus clouds at 78, 000 feet after correction for 
gas absorption; the uncorrected observations are 
inchoated by triangles and dotted bars in the cases 
where they differ significantly from the " corrected" 
values . 

Fig . 9 Same as Fig. 8 for two additional pa rticle size 
distributions . 

Fig. 10 Theoretical cloud reflectivities for 9 = 0°, 9 q = 48 
and ft - ft = 180°, for two cloud optical thicknesses 
for the size distribution of spherical ice parti- 
cles having the characteristic particle size r m = 
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Fig. 11 Theoretical cloud reflectivities for r ^ 10 for five 
size distributions of spherical ice particles. The 
circles represent observations by Blau, et al. of 
" cumulus' 1 clouds at 50,000 feet after correction for 
gas absorption; the uncorrected observations are in- 
dicated by triangles. 

Fig. 12 Theoretical cloud reflectivities as a function of 

azimuth angle for 0 - 80° and 9 = 60° for a size dis- 

o 

tribution of spherical water particles with r m = 2. 

Fig. 13 Same as Fig. 12 for a size distribution with r = 16. 

m . 

Fig. 14 Theoretical and observed cloud reflectivities for 
0 = 80° and 56 ^ 0 Q ^ 64 normalized to unity at 
the smallest observed value of 0 - 0 O * The theoreti- 
cal calculations are for water particles at values 
of 0, 0 Q , and 0 - 0 O , correct for each observed point 
and connected by straight lines. The computations 
were made for t = 32, but are approximately valid for 
t > 10. The observations by Blau, et al. were made on 
cumulus clouds at 4,000 feet above Hurricane Gladys. 

Fig. 15 Theoretical cloud reflectivities as a function of 



azimuth angle for 0 - 80° and 0 O - 50° for a size 
distribution of spherical ice particles having r m = 16|i, 
16 Same as Fig. 15 but for four particle size distributions 
and for only one cloud optical thickness ( t - 32) , 



Fig. 17 Theoretical and observed cloud reflectivities for 
0 = 80° and 44 p ^ 8 Q , £ 56° normalized to unity at 
the smallest observed value of 0 - tf Q . 

Fig. 18 The optical constants for ice are shown in the upper 
part of the figure with the solid curve representing 
the data of Irvine and Pollack (1968 ) r the other two 
curves were obtained by multiplying Irvines and Pollack' 
n r n^ by the factor c in the interval 2.8 ^ \ < 3.6. The 
lower part of the figure shows the single scattering 
parameters < cos 0 > and for the three sets of optical 
constants . 

Fig. 19 Same as Fig. 11 with the theoretical curves for the 
size distribution having r m * 16jj; the three curves 
correspond to the three choices of c shown in Fig. 18. 
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